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Abstract. We discuss a program for systematic studies of heavy quark production in 
pp, pA and AA interactions. The QQ production cross sections themselves cannot be 
accurately predicted to better than 50% at RHIC. For studies of deviations in QQ pro- 
duction such as those by nuclear shadowing and heavy quark energy loss, the pp cross 
section thus needs to be measured. We then show that the ratio of pA to pp dilepton 
mass distributions can provide a measurement of the nuclear gluon distribution. With 
total rates and nuclear shadowing under control it is easier to study energy loss and to 
use cc as a normalization of J/\\t production. 

Keywords: relativistic heavy-ion collisions, heavy flavors 
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1. Introduction 

It is important to have an accurate measure of the charm and bottom cross sections for 
several reasons. Heavy quark decays are expected to dominate the lepton pair continuum 
from the J/\f{cc) and T(bb) up to the mass of the Z° [ 0, §, §]. Thus the Drell-Yan yield 
and any thermal dilepton production will essentially be hidden by the heavy quark decay 
contributions [ |lj|. The shape of the charm and bottom contributions to this continuum 
could be significantly altered by heavy quark energy loss [ |2|, |J. If the loss is large, it may 
be possible to extract a thermal dilepton yield if it cannot be determined by other means 
[ g], Heavy quark production in a quark-gluon plasma has also been predicted [ ||]. This 
additional yield can only be determined if the AA rate can be accurately measured. Finally, 
the total charm rate would be a useful reference for J/\\t production since enhancement of 
the J /\|/ to total charm ratio has been predicted in a number of models [ 0, § [Io[ [llj [j~2|] . 
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2. Baseline Rates in pp 

We first discuss some new calculations of the QQ total cross sections in pp collisions with 
the most recent nucleon parton distribution functions. At leading order (LO) heavy quarks 
are produced by gg fusion and qq annihilation while at next-to-leading order (NLO) qg and 
qg scattering is also included. To any order, the partonic cross section may be expressed in 
terms of dimensionless scaling functions fj-j that depend only on the variable r| [ |T3| ], 

m k=0 1=0 \ m Q ) 

where s is the partonic center of mass energy squared, mg is the heavy quark mass, p is the 
scale and r| = s/Am 2 Q — 1. The cross section is calculated as an expansion in powers of tt s 
with k = corresponding to the Born cross section at order o(oc^). The first correction, 
k=l, corresponds to the NLO cross section at o((Xj). It is only at this order and above 
that the dependence on renormalization scale, pr, enters the calculation since when k = I 
and I — 1, the logarithm ln(/j 2 /mg) appears. The dependence on the factorization scale, 
Pf, the argument of a s , appears already at LO. We assume that pr = pf = p. The next-to- 
next-to-leading order (NNLO) corrections to next-to-next-to-leading logarithm have been 
calculated near threshold [ [l3]] but the complete calculation only exists to NLO. 

The total hadronic cross section is obtained by convoluting the total partonic cross 
section with the parton distribution functions (PDFs) of the initial hadrons, 



{s,m l Q )= £ j^—d(x l x 2 -x)F l l '(x i y)FP(x 2 ,p 2 )a ij (W Q y), (2) 



J- 

'PP 



where the sum i is over all massless partons and x\ and x 2 are fractional momenta. The 
PDFs, denoted by Ff, are evaluated at scale p. All our calculations are fully NLO, apply- 
ing NLO parton distribution functions and the two-loop oc s to both the o(u 2 ) and o((Xj) 
contributions, as is typically done [ [j~3|, p"4[ |. 

To obtain the pp cross sections at RHIC and LHC, we first compare the NLO cross 
sections to the available cc and bb production data by varying the mass, mg, and scale, p, 
to obtain the 'best' agreement with the data for several combinations of mg, p, and PDF. 
We use the recent MRST HO central gluon [ CTEQ 5M [ and GRV 98 HO [ H 
distributions. The results for the cc cross section in pp interactions is shown in Fig. |l . On 
the left-hand side, p = m c for 1.2 < m c < 1.8 GeV, while on the right-hand side, p = 2m c 
for the same masses, all calculated with MRST HO. The scale is not decreased below m c 
because the minimum scale in the PDF is larger than m c /2. The cross sections with p = m c 
are all larger than those with p = 2m c for the same m c because a s (m c ) > a s (2m c ) by virtue 
of the running of a s . Evolution of the PDFs with p tends to go in the opposite direction. At 
higher scales the two effects tend to compensate and reduce the scale dependence but the 
charm quark mass is not large enough for this to occur. 

The best agreement with p = m c is for m c — 1 .4 GeV and m c = 1.2 GeV is the best 
choice for p = 2m c for the MRST HO and CTEQ 5M distributions. The best agreement 
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Fig. 1. Total cc cross sections in pp interactions up to ISR energies as a function of the 
charm quark mass. See [ [l5]] for references to the data. All calculations are fully NLO using 
the MRST HO (central gluon) parton densities. The left-hand plot shows the results with 
p. = m c while in the right-hand plot p = 2m c . From top to bottom the curves are m c = 1 .2, 
1.3, 1.4, 1.5, 1.6, 1.7, and 1.8 GeV. 



with GRV 98 HO is p = m c = 1.3 GeV while the results with p = 2m c lie below the data 
for all m c . All five results agree very well with each other for pp — > cc, as shown on the left 
side of Fig. |[ There is more of a spread in the %~ p ~ + cc results, shown on the right side of 
Fig. g. This is because the K~ PDFs are not very well known. The last evaluations, SMRS 
[ |!|, Owens-7t [ |0|], and GRV-Jt [ |l|] were 10-15 years ago and do not reflect any of the 
latest information on the low x behavior of the proton PDFs, e.g. the distributions are all 
flat as x — > with no low x rise. These pion evaluations also depend on the behavior of the 
proton PDFs used in the original fit, including the value of Aqcd- Thus the pion and proton 
PDFs are generally incompatible. Note that the K~ p cross sections are a bit lower than the 
data compared to the pp cross sections, suggesting that lighter quark masses would tend to 
be favored for this data. The reason is because the low x rise in the proton PDFs depletes 
the gluon density for x > 0.02 relative to a constant at x — ► for p = po the initial scale of 
the PDF. The K~p data are in a relatively large x region, 0.1 < x = 2p/y/s < 0.3, where this 
difference is important. 

We have tried to play the same game with the bb total cross sections but these have 
mostly been measured in K~p interactions. The typical x values of bb production are even 
larger than those for cc but it is not clear that K~p — > bb also favors lower masses. At the 
fixed target energies of bb production, qq annihilation dominates while gg fusion is still 
most important for cc production [ [22]]. The valence-valence u n u p contribution is most 
important since valence distributions dominate at large x. For all three PDFs used, we find 
nij, = p = 4.75 GeV, mj =p/2 = 4.5 GeV, and — 2p — 5 GeV most compatible with the 
sparse data. Attempts to measure the bb total cross section in fixed-target pp interactions 
have been less successful. Hopefully the HERA-B experiment at DESY [ ^3j will soon 
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Fig. 2. Total cc cross sections in pp (left) and n~ p interactions compared to data. See [ |l5|] 
for references (right) to the data. All calculations are fully NLO. The curves are: MRST 
HO (central gluon) with p = m=\A GeV (solid) and fi = 2m = 2.4 GeV (dashed); CTEQ 
5M with n = m=l.4 GeV (dot-dashed) and p = 2m = 2.4 GeV (dotted); and GRV 98 HO 
with /j = m = 1.3 GeV (dot-dot-dot-dashed). 



provide a new measurement. 

Our calculations can then be extrapolated to RHIC and LHC energies. The result for 
cc is shown in Fig. ||. Even though the cross sections agree within 30% at 40 GeV, by 
the Pb+Pb energy of the LHC they differ by a factor of 2.3. The spread in the bb cross 
sections is considerably smaller, ~ 20 — 30% at the ion collider energies. Our results for 
pp interactions at 40 GeV, 200 GeV, and 5.5 TeV are given in Table [j] The AA rates per 
event at b = with the same energies can be obtained by multiplying these cross sections 
by TaaQj = 0), 29.3/mb for Au+Au and 30.4/mb for Pb+Pb. We find 8-13 cc pairs and 
~ 0.05 bb pairs at RHIC with 97-225 cc pairs and ~ 5 bb pairs at LHC without nuclear 
shadowing. The shadowing effect is rather small for cc at RHIC and actually enhances the 
bb rate. The only important modification due to shadowing in the total cross section is on 
the cc rate at the LHC which is reduced to 67-150 pairs. Energy loss does not affect the 
total rate [ Q. As noted by Thews, this cc rate is large enough at RHIC for independent 
c and c quarks to dynamically recombine to form J/\\f's [ pTp . The baseline rates of QQ 
production are thus important for studying these effects. 

3. Nuclear Gluon Distribution in pA 

We now turn to a calculation of the nuclear gluon distribution in pA interactions [ Et|]. We 
show that the dilepton continuum can be used to study nuclear shadowing and reproduces 
the input shadowing function well, in this case, the EKS98 parameterization [ [28} ]. To 
simplify notation, we refer to generic heavy quarks, Q, and heavy-flavored mesons, H. The 
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lepton pair production cross section is 

dM dy = J d3 Pl d3 PlJ d3 PHd 3 PH^M n -M(p,,p I ))8(y ll -y(p hP j)) 

dT H ^ l+x {p H ) dTB^+Xjpjj) d^ m+x 
(Ppi cPpj cfipucPpfj 

X6(ymjn < yijj < ymax)8(<l>iiiin < < <l>max) (3) 

where M(pi,pj) and y(pi,pj) are the invariant mass and rapidity of the 11 pair. The decay 
rate, dT H ^ l+x (pn)/(Ppi, is the probability that meson H with momentum pn decays to 
a lepton I with momentum pi. The 9 functions define single lepton rapidity and azimuthal 
angle cuts used to simulate detector acceptances. 

Using a fragmentation function Dq to describe quark fragmentation to mesons, the 
HH production cross section can be written as 

da pA^HH+X t&fn&f7i doP A ^QQ+ X f H N 77, 

x ^XpH-ZiPQ^iPTj-ZiPQ). (4) 

Our calculations were done with two different fragmentation functions. We found that our 
results were independent of Dq. The hadronic heavy quark production cross section per 
nucleon in pA collisions can be factorized into the general form 

1 dc pA ^ Q ® +x f d6 ij ^ Q Q 
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Table 1. Charm and bottom total cross sections per nucleon for the extrapolated calcula- 
tions shown previously. The heavy quark mass and factorization/renormalization scales are 
given, along with the cross sections at 40 GeV (HERA-B), 200 GeV (Au+Au at RHIC), 
and 5.5 TeV (Pb+Pb at LHC). 









40 GeV 


200 GeV 


5.5 TeV 


PDF 


m c (GeV) 


cc 
n/m c 


a Cub) 


a C"b) 


a (mb) 


MRST HO 


1.4 


1 


37.8 


298 


3.18 


MRST HO 


1.2 


2 


43.0 


382 


5.83 


CTEQ 5M 


1.4 


1 


40.3 


366 


4.52 


CTEQ 5M 


1.2 


2 


44.5 


445 


7.39 


GRV 98 HO 


1.3 


1 


34.9 


289 


4.59 






bb 








PDF 


m b (GeV) 


fi/mb 


a(nb) 


oC"b) 


oC"b) 


MRST HO 


4.75 


1 


9.82 


1.90 


185.2 


MRST HO 


4.5 


2 


8.73 


1.72 


193.2 


MRST HO 


5.0 


0.5 


10.96 


2.16 


184.8 


GRV 98 HO 


4.75 


1 


13.40 


1.65 


177.6 


GRV 98 HO 


4.5 


2 


12.10 


1.64 


199.0 


GRV 98 HO 


5.0 


0.5 


14.80 


1.73 


166.0 



where ff = Ff jx and ff — Ff jx with Ff — FfRf. The shadowing ratio Rf is that of 
EKS98 [ ^8||. The partonic cross section is the differential of Eq.jjl]) at k = 0. Note that the 
total lepton pair production cross section is equal to the total QQ cross section multiplied 
by the square of the lepton branching ratio. 

We compare the ratios of lepton pair cross sections with the input R^ in Fig. ||. All 
the results are integrated over the rapidity intervals appropriate to the PHENIX and ALICE 
dilepton coverages. The ratio follows R^ at all energies. The higher the energy, the better 
the agreement: at the LHC the two agree very well. 

The ratio always lies below R^ for two reasons. First, qq annihilation is included and 
quark shadowing is different than gluon shadowing. The qq contribution decreases with 
energy, leading to better agreement at the LHC. Second, the phase space integration smears 
the shadowing effect relative to Rg((x2), (fi)). Note that the ratio deviates slightly more 
from R™({x2), (p)) for e + e~ than for fi + ff~ because the curvature of Rg with x is stronger 
at larger values of x and, due to the differences in rapidity coverage, the average values of 
X2 are larger for e + e~ . 

The average X2 decreases with energy. We have 0.14 < (X2) < 0.32 at the SPS where 
Rg is decreasing. At RHIC, 0.003 < (X2) < 0.012, where R^ is increasing quite rapidly. 
Finally, at the LHC, 3 x 10~ 5 < (^2) < 2 x 10~ 4 where is almost independent of x. The 
values of (^2) are typically larger for electron pairs at collider energies because the electron 
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Fig. 4. The ratios of lepton pairs from correlated DD and BB decays in pA to pp colli- 
sions at the same energies (solid curves) compared to the input R% at the average xi and p 
(dashed)/ y^u 5 ) (dot-dashed) of each M bin. From Ref. [ 0. 

coverage is more central than the muon coverage. 

The average p? increases with energy and quark mass. For cc we have 7.58 < (p 2 ) < 
48.5 GeV 2 at the SPS, 9.46 < (p 2 ) < 141 GeV 2 at RHIC, and 11.4 < (p 2 ) < 577 GeV 2 
at the LHC. For bb production, 32.0 < (p 2 ) < 54.3 GeV 2 at RHIC and 37.9 < {p 2 ) < 156 
GeV 2 at LHC. 

4. Heavy Quarks in AA 

4.1. Effects of Energy Loss 

Energy loss would best be determined by reconstruction of D and B meson decays and 
comparing with distributions expected from pp and pA extrapolations that do not consider 
energy loss. Whether energy loss is measurable in reconstructed D and B decays or not, the 
change in the dilepton continuum should surely be present if the loss is nonzero and will 
bias the interpretation of the dilepton continuum. So far, the amount of the energy lost by 
heavy quarks is unknown. While a number of calculations have been made of the collisional 
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loss in a quark-gluon plasma [ g^], only recently has radiative loss been applied to heavy 
quarks [ |3^|. The radiative loss can be rather large, dE/dx ~ —5 GeV/fm for a 10 GeV 
heavy quark, and increasing with energy, but the collisional loss is smaller, dE/dx ~ 1 — 2 
GeV/fm, and nearly independent of energy [30]. We note that energy loss will suppress 
high pr and large invariant mass quark pairs as long as \dE/dx\ > (pr) /Ra [ [ji- 
lt is important to note that energy loss does not reduce the number of QQ pairs pro- 
duced but only changes their momentum. However, an effective reduction in the observed 
heavy quark yield can be expected in a finite acceptance detector because fewer leptons 
from the subsequent decays of the heavy quarks will pass kinematic cuts such as a mini- 
mum lepton pr. 

If the loss or the pj cut is large, the Drell-Yan and thermal dileptons could emerge 
from under the reduced DD and BB decay contributions at large masses. Even without con- 
sidering energy loss, Gallmeister et al. suggested that thermal dileptons could be detected 
by increasing the minimum lepton pj because, in the D and B rest frames, the maximum 
energy of the individual leptons is limited to 0.9 and 2.2 GeV respectively. The lepton pj 
from thermal production has no such limitation [ B] . 



4.2. Quarkonium normalization 

Heavy quark production in AA collisions is also interesting because of the prominent ef- 
fect it could have on quarkonium. Initial nucleon-nucleon collisions may not be the only 
source of quarkonium production. Regeneration of quarkonium in the plasma phase [ 
0, ||, H, [Hj [jr], Q could counter the effects of suppression, ultimately leading to en- 
hanced quarkonium production. In the plasma phase, there are two basic approaches: sta- 
tistical and dynamical coalescence. Both these approaches depend on being able to measure 
the quarkonium rate relative to total QQ production. Thus the QQ rate is preferable as a 
normalization of quarkonium production, particularly since both share the same production 
mechanisms and approximate (x), (jj) values. However the final-state effects such as energy 
loss will make the total rate difficult to quantify without substantial detailed studies. These 
secondary production models should be testable already at RHIC where enhancements of 
factors of 2-3 are expected from coalescence [ ^ [Tl| ]. 

Other processes besides heavy quark production have been suggested as references for 
quarkonium production. Using the Z° as a reference [ |3l|] as a reference would eliminate the 
uncertainty due to final-state effects on the Z° — > l + l~ decays but the different production 
mechanisms and masses leaves it less desirable. It has also been suggested that the \(/'///\|/ 
and T'/T ratios be studied as a function of pj [ |32|] since deviations from the pp ratios 
should reflect quark-gluon plasma characteristics. The only drawback to such a mechanism 
is that strong suppression may result in poor statistics. 
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